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Abstract. In this paper, we introduce 3-dimensional L— summing method, which is a rearrangement of 
the summation "^Aatc with I < a,b, c < n. Applying this method on some special arrays, we obtain some 
identities on the Riemann zeta function and digamma function. Also, we give a Maple program for this 
method to obtain identities with input various arrays and out put identities concerning some elementary 
functions and hypergeometric functions. Finally, we introduce a further generalization of L— summing 
method in higher dimension spaces. 



1. Introduction and Motivation 
Consider the following n x n multiplication table 



1 


2 




k 




ri 




k 


2 


4 




2k 




2)7 


2k 
















A- 


2k 




k-' 




ku 


k 2k • • ■ 


















>i 


•2n 




kn 




■> 
»- 



Figure 1. Multiplication table and L— summing element, Lk 



If we set S(7T.) for the sum of all numbers in it, then by summing line by line we have — ^ "("^+^) j ^ Qn 
the other hand, we can find S(n) by using another method; letting be the sum of numbers in the rotated 
L in above table (right part of FIGURE 1), we have 

Lfe = fc + 2fc H h H \-2k + k = 2fc(l + 24 h fc) - fc^ = fc^. 

We call Lk, L— summing clement. Thus we get S(7i) = J2k=i^k = J^k^i^^' ^^"^ therefore X]fc=i — 
^ n(rt^+i) ^ This is 2-dimensional L-summing method (applied on the array Aat = ab), which briefly is 

(1.1) Z^i^ ~ Summing Elements) = E. 
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More precisely, the L— summing method of elements of n x n array Aab with 1 < a,b < n, is the following 
rearrangement 



fc 




AkA= E ^ 



ab • 

l<a.b<n 

This method allows to obtain easily some classical algebraic identities and also, with help of MAPLE, 
some new compact formulas for sums related with the Riemann zeta function, the gamma function and the 
digamma function [SJ [3] . 

In this paper we introduce a 3-dimensional version of L— summing method for nx n x n arrays and applying 
it on some special arrays we obtain some identities concerning the Riemann zeta function and digamma 
function. Then, we give a Maple program for this method and using it we generate and then proof some 
new identities, concerning some elementary functions and hypergeometric functions. Finally, we introduce 
a further generalization of L— summing method in higher dimension spaces and for latices related by a 
manifold. 

2. Formulation of the L— summing method in 
Consider a three dimensional array Aabc with I < a,b,c < n and n is a positive integer. We should 



prepare an explicit version of the general formulation (1.1) for this array. The summation of all entries is 



S(") = X^Ka b c<n Aabc- Thc L— Summing elements in this array have the form pictured bellow 




Figure 2. L— summine; elements in 



So, we have = E2 — Si + Sq, with 

k k k k k k 

S2 = E Aki,c.+ E Aakc+ E Aabk, ^1 —''^Aakk + ''^Akbk + ''^Akkc, ^ Akkk- 

6,c— 1 a.c—1 a, 6—1 a—1 b—1 c—1 

Therefore, L— summing method in M'^ take the following formulation 

n 

(2.1) E{S2-Si+Eo} = S(n). 

fc=i 

Note that S2 is the sum of entries in three faces, Ei is the sum of entries in three intersected edges and So 
is the end point of all faces and edges. 
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If the array Aabc is symmetric, that is for each permutation a Cz S3 it satisfies Aabc — ^o-(a)CT(fc)CT(c)i then 
L— summing elements in takes the following easier form 



k 

(2.2) ife = 3 ^ Akbc ~ 3 Agkk + ^fcfcfc- 

l<6,c<fe a=l 

In the next section we will apply 3-dimensional L— summing method on two special symmetric arrays, related 
by the Riemann zeta function and digamma function. 

3. Arrays related by the Riemann zeta function and digamma function 

3.1. The Riemann zeta function. Suppose s G C and let Aabc — {abc)^'^. It is clear that 

3 



l<a,b,c<n \k=l / 



where Cn{s) — X]fc=i ^ Since this array is symmetric, considering (2.2 1, we have 
and an easy simplifying, we can reform J2^k — as follows 

(^-,. ^CU^) Ckjs) ^ CUs)-U3s) 

fc=i 

Note that if 5R(s) > 1, then lim„^oc Cn{s) — C{s), where C(s) = X^fcLi is the well-known Riemann zeta 
function defined for complex values of s with 3fi(s) > 1 and admits a meromorphic continuation to whole 
complex plan [5]. So, for 5ft(s) > 1 we have 

^ fc^ fc2s 3 

fe=l 

which also is true for other values of s by meromorphic continuation, except s — 1 and s = ^. 



3.2. Digamma function. Setting s = 1 in (3.11 (or equivalently taking Aabc — ^) and considering 

1 

fe' 

E 



Cn(l) = = Ylk=i h obtain 



Hi H, _ Hi - C„(3) 



k k"^ 3 ' 

One can state this identity in sense of digamma function ^{x) = ^ lnr(a;), with r(a;) = e~^t''^~^dt is 
the well-known gamma function. Considering logarithmic derivative of the formula r(n + 1) = nT(n), we 
obtain 

(3.2) *(n-M) = - + *(n), 

and applying this relation, we yield that ^'(n + 1) — ^'(1) = Hn- Thus, we have 

(3.3) ^'(n+l) + 7 = i/„. 
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in which 7 = 0.57721 • • • is the Eulcr constant [1 . Therefore, we obtain 



(3.4) 

Letting 



E 

fe=i 



(^(fc+ 1) +7)^ _ ^(fc+ 1) + 7 ) _ {^{n + 1) + 7)^ - Cn(3) 
k 



6{m,n) = ^ 



k=l 



the following identity ^ is a result of 2-diniensional L— summing method 
(3.5) 



e(l, n) = - - - + 1)^ - 



where ^(m, x) = ^^-^^(x) is called m*'' polygamma function [T and we have X]fc=i w = ^ ^(1, + 1) + if 
which is special case of the following identity 



(3.6) 



C«(s)-|^^*(.5-l,n+l) + C(s) 



(s e Z,s > 2), 



and using it in (3.1 1 one can get a generalization of (3.4), however this relation itself is the key of getting an 
analogue of (3.5 1 in K"^, stated bellow. 

Theorem 3.1. For every integer n > I, we have 



E 

k=l 



*(fc)2 xSj{k) (*(n+l) + 7)3 C«(3) 



fc2 



+ (7-2)— -(7-2)vI/(l,n+l)-72vI/(n + l)-73-26(l,n). 
3 6 



Proof. We begin from the left hand side of the identity (3.4 1; simplifying it by the relations (3.2 1, (3.3 1 and 



the relation (3.6) with s = 2, gives the result. 

Corollary 3.2. For every integer n > I, we have 
(vI/(n+l) + 7)3 C«(3) 



□ 



6(2, n) 



(7-2)y- (7-2)^(1,71+1) -72vl/(n+l) -73 -26(1, n)-^ 



^ *(fc) 



fe=i 



fc2 



In above corollary, the main term in the right hand side is 



■i/{n+l) 



Also, we note that the summation 



X]fc=i ^^fef^ converges. Thus, we can write the following asymptotic relation 



6(2, n) = 



^'(n+ 1)' 



+ 0{ln^ n) 



Similarly, considering (3.5) we have 



^(n + 1^)2 
6(1, n) = ^ ^ ' + 0(ln n) 



{n —^ 00). 



(n — > 00). 



Note and Problem. It is interesting to find an explicit (probably recurrence) relation for the function 6(m, n). 
Considering two above asymptotic relations, we guess that 

&(m,n) — \- Oiln™ n) (n— s-cxd). 

^ ^ m + 1 ^ ^ ^ ^ 

One can attack to this problem considering generalization of L— summing method in higher dimension spaces, 

pointed in the last section of this paper. 
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4. Generating some new identities by Maple and 3-dimension i— summing method 

Appendix of this paper includes Maple program of 3-dimension L— summing method. By LSMI < Aatc >, 
we call the identity outputted by L— Summing method's Maple program with input Aatc- The algorithm 
of this program is result of the formulation of 3-dimension L— summing method in above sections. In this 
program we input a 3-dimensional array Aabc, then out put is an identity generated by Maple. In this section 
we will state some of these identities, with handling a detailed proof. 

4.1. Some elementary functions. 

Proposition 4.1. We have 

n 

LSMI<ln{a) >: ^ {/c^ Infc + 2 fclnr(fc + 1) - 2 fcln A: - lnr(fc + 1) + In A:} = n^lnr(n+ 1). 
fc=i 



Proof. Considering the relations (2.11 and T{n + 1) = nl, we have T,{n) = "-^X]a=i^^'^ ~ lnT{n + 1). 
Also, E2 = fc2lnfc + 2fclnr(fc-|- 1), Ei = lnr(fc + 1) + 2klnk and Eq = Infc. Putting these relations in ( |2l" l 
yields LSMI < ln(a) > as desired. □ 

Corollary 4.2. We have 

n 

^{(fc2-/c)lnfc-K2/clnr(fc+l)} ^ (n^ + n) lnr(n + 1). 
fc=i 



Proof. Breaking up the statement under the summation obtained from LSMI < ln(a) > in Proposition 4.1 
into the sum of (fc^ — /c) In fc + 2fc In r(fc + 1) and In F (fc + 1) + In fc — In fc, and considering the Proposition 
6 of [3], which states 

n 

^{lnr(A: + 1) + fclnfc - Infc} = nlnr(n+ 1), 
fc=i 

completes the proof. □ 

Remark 4.3. Examining Maple code of expressed summation on above corollary, one can see that Maple 
has no comment on the computing this summation; however, it is obtained by Maple itself and L— summing 
method. This example shows that program-writers of Maple can add L— summing method in the summation 
package of this software, in order to making it able to compute some summations which already couldn't 
compute them. 

Proposition 4.4. A little simplifying LSMI < tan(a) >, we have 

n 

{(fc- l)2tanfc+ (2fc - l)T(fc)} = nH{n), 

k=l 

where T(n) — J2k=i ^• 



Proof. Considering the relation (2.1 1, we have S(n) — X]a=i tana = n'^1{n). Also, S2 = fc^ tan fc+2fcT(fc), 
El = T(fc) + 2fctanfc and Eq = tanfc, and consequently Lk — {k — l)^tanfc + (2fc — l)X(fc). This completes 
the proof. □ 
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4.2. Hypergeometric functions. In the next proposition, we introduce an identity concerning hyperge- 
ometric functions, denoted in Maple by hypergeomC [oi 02 • • • ap] , [&i 62 • • • 6g] ,x). Standard notation 
and definition [B] is as follows 



ai 02 



k>0 



where 



+ 1 



bi 62 • • • bg 

{k + ai){k + 02) • • • (fc + flp) 



tk {k + bi){k + b2)---{k + bq){k + l) 
Proposition 4.5. A little simplifying LSMI < a! > and stating in standard notations, we have 

n 

{(fc- l)2fc! + (2fc- l)(fc+ l!)i3(l,fc + 2)} = l)!i3(l,n + 2), 



fe=i 



whe 



a P 



Proof. Considering definition of hypergeometric functions we have n + I) — (n + 1)^)(1, n + 2), which 
implies J22=i ~ 2) — (^^ + n + 2) = say. This gives S(7i) = and in similar way 

it yields that Lk = (fc - l)^fc! + (2fc - 1) ((fc + fc + 2) - ^(1,2)). Thus, we obtain 

n 

^ {(fc - l)^fc! + (2fc - 1) ((fc + l\)9j{l, fc + 2) - fi(l, 2))} = n^{n + iy.Sj{l,n + 2) - ^^^^(l, 2), 



fe=i 



and a easy simplifying this, implies the result. 



□ 



Remark 4.6. Three last propositions are examples of the array Aabc = f{oL), for some given function /. In 
this case, L— summing method takes the following form 



k=l 



where ^{n) = YZ=i /(«)■ 



5. Further generalizations of the L-summing method and some comments 

5.1. The i— Summing method in M*. Consider a i— dimensional array Ax^x2---xt ^^nd let S(n) — ^ Ax^x2---xt 
with 1 < xi,X2,--- ,xt < n. The i— Summing method in M* is the rearrangement S(n.) = ^L^, where 
Lk = EL=i {(-1)"-^S,_„} and 



I<il<i2<---<jm<i 



where in the inner summation is over Xj g {xij,-- - ,2;^^}*^ = {a;i,a;2j''' j^j^} — {xij^,--- ,Xi^'\ with 
1 < < fc, and the index Xiji^ - im denotes a;ia;2 ■ ■ ■ Xt with a;ij^ = = ■ ■ ■ — — fc. One can apply 
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this generalized version to get more general form of relations obtained in previous sections. For example, 
considering the array A^-^^^...^^ — {xiX2 ■ • • Xt)^^ with s G C, yields 

E I E (-1)""' Cm) ^"™'^^(^)*"™| = ^"(■'*)* + i-^yuts). 

k=l lm=l ^"^^ J 

5.2. L— summing method on manifolds. As in rising of this paper, the base of the i— summing method is 
ordinary multiplications table. Above generalization of the L— Summing method in M* is based on generalized 
multiplication tables [1]. But, is a very special <-dimensional manifold, and if we replace it with F, an 
I— dimensional manifold with / < t, then we can define generalized multiplication table on F by considering 
lattice points on it (which of course isn't easy problem). Let 

Lr{n) = {(ai, a2, • • ■ , at) e F n N* : 1 < ai, a2, • • • ,at < n}, 

and / : M'^' — > C is a function. If Or is a collection of /c — 1 dimension orthogonal manifolds, in which 
Ly^ti) = UAeOr-^A(^) and Ljy.(n)nL/^^(n) — (j) for distinct A^, Aj G Op, then we can formulate L— summing 
method as follows, 

E /w = E ( E /(^) 

XeLr{n) AeOr [xsLaC") 

Here L— summing elements are X]jfeLA(n) /("^)- This may be useful when one apply it on some special 
manifolds. 

5.3. Stronger form of L— summing method. One can state the method of i— summing 'Y^L^ = 2](n) 
in the following stronger form 

Ln = S(n) — S(n — 1). 

Specially, this will be useful for those arrays with S(n) computable explicitly and Lk maybe note. For 
example, applying this note on the array Ax-^x2---xt — {xiX2 ■ ■ ■ Xt)^^ in M* with s G C, implies 



(-l)'"-'Q)n-'-Cn(s)*-™ = UsY + {-l)\n{ts) - Cn-l{sY - {-if Cn-l{ts) . 



t-1 

E^ -"-^ 

m— 1 
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Appendix. Maple Program of 3-dimension i— Summing Method for the array Aabc — ^ 
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restart : 

A [abc] :=l/(a*b*c) ; 

S21 :=siim(svim(eval(A[abc] ,a=k) ,b=l. .k) ,c=l. .k) : 

S22:=sum(sum(eval(A[abc] ,b=k) ,a=l. .k) ,c=l. .k) : 
S23:=sum(sum(eval(A[abc] ,c=k) ,a=l. .k) ,b=l. .k) : 
S2:=S21+S22+S23: 

Sll:=siiin(eval(eval(A[abc] ,a=k) ,b=k) ,c=l. .k) : 
S12:=sum(eval(eval(A[abc] ,a=k) ,c=k) ,b=l. .k) : 
S13 : =sum(eval(eval(A [abc] ,b=k) ,c=k) ,a=l . .k) : 
S1:=S11+S12+S13: 

SO : =eval (eval (eval (A [abc] ,a=k) ,b=k) ,c=k) : 
L[k] —simplify (S2-S1+S0) : 

ST (A) := (simplify (sum (sum (sum (A [abc] ,a=l . .n) ,b=l . .n) ,c=l . .n) ) ) : 
Sum(L[k] ,k=l. .n)=ST(A) ; 

A, 

abc 

(4'(fc + l))"fe^ + 6*(fc + l)fc^7 + 372/c2-3*(/c + l)fe-37fe + l , ,3 
p = (* (n + 1) + 7) 

fe=i 
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